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1.  Introduction 


A  potential  consequence  of  a  Titan  IV  failure  early  in  flight  is  the  ground  impact  of  a  Solid 
Rocket  Motor  Upgrade  (SRMU)  segment.  Under  certain  meteorological  conditions,  the 
explosion  of  propellant  which  results  from  impact  may  generate  high  overpressures  (enough 
to  shatter  windows)  many  miles  away  from  the  blast  source.  BLASTC  (Ref.  1)  is  a  tool 
developed  by  ACTA  Inc.  for  predicting  window  breakage  and  casualties  resulting  from  such 
launch  failure  overpressures  at  Vandenberg  Air  Force  Base. 

The  yield  histogram  input  to  BLASTC  currently  assumes  that  all  SRMU  segment  impacts  are 
side-on.  The  rationale  behind  this  assumption  is  as  follows.  Hydro-code  calculations  of  yield, 
generated  for  perfectly  side-on  and  perfectly  end-on  impacts  at  the  same  velocity,  show  that 
the  side-on  case  gives  significantly  higher  yield  (Ref.  2.)  There  are  no  code  results  (nor  test 
data)  at  intermediate  impact  angles.  Furthermore,  it  had  been  previously  suggested  to  ACTA 
by  another  agency  that  free-falling  cylindrical  bodies  would  tend  to  aerodynamically  trim  to 
side-on  orientation  prior  to  impact.  Thus,  lacking  further  information,  ACTA  adopted  the 
logically  conservative  approach  of  assuming  side-on  impacts  in  all  cases. 

Since  the  assumption  of  side-on  impact  for  all  failures  may  be  a  source  of  undue  conservatism, 
the  following  two  tasks  were  undertaken  by  the  Aerospace  Flight  Mechanics  Department. 

1.  Calculation  of  a  probability  distribution  of  impact  orientation  angles  using  Monte  Carlo 
techniques  to  randomly  sample  initial  conditions  at  some  randomly  selected  time  of 
failure.  For  given  initial  conditions,  the  free-falling  six  degree-of-freedom  trajectory 
is  simulated  to  the  ground,  using  realistic  aerodynamic  forces  and  moments  on  the 
SRMU  segments.  The  impact  angle  in  each  case  is  recorded. 

2.  Formulation  of  a  method  by  which  a  given  impact  orientation  angle  can  be  associated 
with  a  yield  value  lying  in  the  range  between  the  side-on  and  end-on  values. 
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2.  Probability  Distribution  of  Impact  Angles 


A  Monte  Carlo  trajectory  simulation  technique  was  used  for  approximating  the  segment 
impact  orientation  distribution.  The  simulation  initiates  a  nominal  Titan  IVB-12  three 
degree-of-freedom  (3D0F)  performance  trajectory.  At  a  selected  time  of  failure,  the  simula¬ 
tion  “jumps”  to  modeling  an  SRMU  segment  in  a  6DOF  free-fall  (non-thrusting)  to  ground 
impact.  The  mass  and  moments  of  inertia  are  calculated  based  on  the  elapsed  trajectory 
time.  Aerodynamic  coefficients  of  cylindrical  segments  of  various  lengths  were  provided  by 
the  Aerospace  Fluid  Mechanics  Department  (Ref.  3.)  These  coefficients  are  shown  as  func¬ 
tions  of  angle-of-attack  in  Figs.  1  -  12.  Discussions  with  ACTA  (Ref.  4)  revealed  that 
scenarios  which  result  in  an  intact  segment  impact  assume  a  failure  (or  partial  failure)  of 
the  Inadvertent  Separation  Destruct  System  (ISDS.)  Such  a  condition  is  associated  with  an 
SRMU  case  rupture  on  an  otherwise  nominal  flight.  This  would  imply  an  on-trajectory  failure 
with  initial  angle-of-attack  and  angular  rates  near  nominal.  Given  these  factors,  the  Monte 
Carlo  analysis  assumed  a  uniformly  random  failure  time  in  the  first  60  seconds  of  flight,  and 
random  normal  variation  of  angle-of-attack  and  pitch/yaw  rates  about  the  nominal  state. 
A  standard  deviation  of  5°  was  chosen  for  the  angle-of-attack  probability  distribution.  The 
standard  deviation  for  the  pitch  and  yaw  rates  were  both  chosen  to  be:  30  deg/sec  for  a 
high  rate  analysis,  10  deg/sec  for  a  medium  rate  analysis,  and  3  deg/sec  for  a  low  rate  analy¬ 
sis.  Another  consideration  is  the  possibility  of  two  segments  still  joined  together  at  impact. 
Thus,  cylinder  lengths  of  16,  33,  48,  and  68  ft  were  each  modeled  separately.  (In  all  cases, 
the  diameter  is  10.5  ft.) 

There  were  therefore  a  total  of  twelve  cases  considered,  viz.,  the  high,  medium,  and  low 
angular  rates  for  each  of  the  four  lengths.  In  each  of  these  cases,  the  Monte  Carlo  sample 
size  was  1000.  For  each  element  in  a  sample,  the  failure  time  was  drawn  uniformly  from  the 
first  60  seconds  of  flight.  Variations  about  the  nominal  trajectory  values  (at  time  of  failure) 
of  angle-of-attack  and  rotation  rates  were  realized  by  drawing  from  normal  distributions 
having  the  standard  deviations  described  above. 

The  Monte  Carlo  results  for  the  twelve  cases  are  shown  in  Table  1,  along  with  the  corre¬ 
sponding  impact  orientation  probabilities.  Impact  orientation  angle  histograms  for  the  high 
and  low  rate  analyses  for  the  various  segment  lengths  are  shown  in  Figs.  13-20.  [Histograms 
for  the  medium  rate  (10  deg/sec)  analyses  are  not  shown,  but  they  look  much  like  much  like 
those  for  the  low  rate  (3  deg/sec).] 


Observations 

1.  From  the  histograms  for  the  four  cases  with  the  higher  standard  deviations  in  initial 
angular  velocity  (30  deg/sec),  it  can  be  seen  that  the  orientation  at  impact  tends  to 
be  uniformly  distributed.  This  results  from  the  fact  that  for  higher  initial  angular 
velocities,  there  is  insufficient  time  in  the  cases  considered  here  for  the  body  to  trim 
due  to  aerodynamic  moments. 
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2.  For  lower  initial  angular  rates  ( apy  of  10  deg/sec  and  3  deg/sec),  the  longer  cylinders 
(two  joined  segments,  either  48  or  68  ft  in  length)  trim  to  side-on  attitude  at  impact 
in  slight  preference  to  end-on  attitude. 

3.  For  lower  initial  rates,  the  short  segment  (16  ft)  tends  to  trim  to  an  end-on  impact, 
while  the  single  segment  (33  ft)  tends  to  prefer  a  nearly  end-on  impact  (70°  to  80°.) 
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c Tpy  =  one  sigma  value  for  initial  pitch/yaw  body  rates  (deg/sec) 

9  =  impact  orientation  angle  from  horizontal  (deg) 
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3.  Association  of  Impact  Angle  with  Yield  Value 


As  stated  previously,  there  are  no  code  results  or  test  data  for  the  explosive  yield  resulting 
from  the  impact  of  propellant  segments  at  orientations  other  than  perfectly  end-on  or  per¬ 
fectly  side-on.  In  what  follows,  a  method  is  proposed  by  which  a  yield  value  may  be  assigned 
to  an  impact  of  arbitrary  orientation. 


Mechanical  Rotation  Time 

Consider  a  cylindrical  propellant  segment  of  length  L  and  diameter  D  impacting  the  ground. 
The  impact  angle,  i.e.,  the  angle  between  the  cylinder  axis  and  the  ground  plane,  is  denoted 
by  9  (see  Fig.  21.)  It  is  assumed  that  the  point  of  impact  (shown  as  point  O  in  Fig.  21) 
instantaneously  becomes  a  center  of  rotation,  i.e.,  the  segment  is  assumed  not  to  bounce  on 
impact.  Let  v\  denote  the  vertical  component  of  center-of-mass  velocity  just  prior  to  impact, 
and  let  v2  denote  the  value  of  this  component  just  after  impact. 

Conservation  of  the  vertical  component  of  linear  momentum  implies 

Jl  =  (vi-v2)m  (1) 

where  Jl  is  the  linear  impulse  delivered  by  the  ground  force  during  impact  and  m  denotes 
the  mass  of  the  segment. 

The  angular  impulse  due  to  the  force  from  the  ground  induces  a  post-impact  angular  velocity, 
u>2,  which  is  assumed  here  to  be  very  large  in  comparison  to  any  pre- impact  angular  velocity. 
Conservation  of  angular  momentum  implies 

Ja  —  hu 2  (2) 

where  Ja  is  the  angular  momentum  of  the  ground  force  about  the  center  of  mass  and 

mL 2 

c  _  12 

is  the  moment  of  inertia  of  the  cylinder  about  a  central  diameter. 

From  the  rotational  kinematics  of  the  problem, 


II 

5* 

(4) 

and 

V2  —  lx 

(5) 

where 

lx  =  l  cos  6' 

(6) 

and,  from  Fig.  21 

l  =  \  VL2  +  D2  and  6'  =  0  +  tan ~l(D/L) 

(7) 

1  + 


(3) 
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Substitution  of  (1),  (4),  and  (5)  into  (2)  eliminates  impulse  and  gives  a  single  equation  for 
0>2. 

lx  ^1 


U>2  = 


(8) 


Ic/m  + 1\ 

Now,  the  time  for  the  cylinder  to  rotate  into  full  lateral-tangent  contact  with  the  ground  is 


6 

At  rot  “ 

00  2 


(9) 


Combustion  Wave  Time  Scale 

It  is  now  desired  to  compare  the  rotation  time  in  (9)  to  the  time  taken  for  a  combustion 
wave,  initiated  by  contact  at  point  O  (see  Fig.  21),  to  propagate  the  axial  extent  L  of  the 
segment. 

It  is  assumed,  simplistically,  that  the  combustion  wave  initiated  by  ground  contact  propa¬ 
gates  vertically  to  the  top  surface,  then  reflects  in  accord  with  Shell’s  law.  This  process  is 
repeated  each  time  the  wave  front  reaches  the  lateral  surface  of  the  cylinder.  Thus,  the  axial 
component  of  the  propagation  velocity  of  the  wave,  due  to  multiple  internal  reflections,  is 

cs  sin# 

where  cs  is  the  combustion  wave  speed  in  the  propellant.  The  combustion  time  scale  is 
therefore 

Atc  =  -tA- 
cs  sm# 

and  side-on  impacts  may  now  be  defined  as  those  for  which 

A  trot  <  Atc  (11) 

That  is,  side-on  impacts  are  defined  as  those  cases  in  which  the  cylinder  rotates  into  full 
lateral-tangent  contact  with  the  ground  prior  to  the  creation  of  an  axial  relief  wave  caused 
by  incident  wave  reflection  from  the  end  of  the  cylinder  most  distant  from  the  initiation. 

Substitution  of  (8),  (9),  and  (10)  into  (11)  gives 

( #  sin#)  ( +  lx  )  <1  (12) 

v\L  \  lx  ) 

as  the  condition  to  be  met  in  order  for  an  impact  to  be  considered  side-on. 

The  deflagration  wave  speed,  cs,  in  the  propellant  may  be  estimated  with  the  aid  of  data 
which  appear  in  Fig.  4  on  page  23  of  Ref.  2.  These  data  are  from  the  SOPHY  test  series  and 
were  taken  on  propellant  segments  which  were  six  feet  in  diameter.  Other  data  are  shown  for 
smaller  diameter  segments,  but  the  SRMU  is  10.5  feet  in  diameter,  which  makes  the  SOPHY 
data  the  most  pertinent.  As  stated  on  page  22  of  Ref.  2,  the  composition  of  the  SOPHY 
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propellant  is  similar  enough  to  the  SRMU  propellant  that  a  single  reaction  model  suffices  for 
both.  Averaging  the  13  measured  deflagration  velocities  which  were  taken  at  propagation 
distances  between  1.5  and  4.0  diameters  from  the  booster  charge  gives 

cs  =  3.17  km/sec  =  10,400  ft/sec 

as  the  average  deflagration  speed.  (By  way  of  comparison,  the  Chapman- Jouget  detonation 
speed  in  TNT  is  4.2  km/sec.)  With  cs  so  determined,  equation  (12)  is  now  used  to  plot 
the  threshold  angle,  9C,  below  which  an  impact  may  be  considered  to  be  side-on.  (See  Figs. 
22-25.) 


Assignment  of  Yield  to  an  Arbitrary  Impact  Orientation 

It  is  now  desired  to  assign  a  yield  value  to  an  impact  for  which  the  speed,  V\,  and  angle,  9, 
are  given.  It  is  assumed  that  values  for  the  side-on  yield,  zSide,  and  the  end-on  yield,  zend->  are 
known.  The  following  procedure  is  proposed.  For  the  segment  length  and  impact  velocity 
at  hand,  the  threshold  angle,  0C,  is  found  from  the  plots  in  Figs.  22-25.  The  yield  is  then 
assigned  to  one  of  the  following  three  values. 


'  Zside  if  9  <  9C 

<  Zend  if  9  >  2  9C 

( Zend  +  Zside)/ 2  otherwise 


The  physical  reasons  for  this  assignment  are,  that  in  the  first  case,  the  segment  rotates  into 
full  lateral  contact  with  the  ground  before  the  axial  deflagration  can  propagate  the  length 
of  the  segment,  implying  that  excitation  is  dominated  by  side-on  impact.  In  the  second 
case,  the  axial  deflagration  has  had  the  opportunity  to  propagate  the  length  of  the  segment 
at  least  twice  before  full  side  contact  is  made  with  the  ground.  Subsequent  excitation  by 
side-on  contact  is  secondary,  and  the  explosion  strength  will  be  more  directly  dictated  by 
ground  contact  with  the  end  of  the  cylinder.  The  third  case  is  the  intermediate  situation  in 
which  a  simple  average  of  the  end-on  and  side-on  values  is  used. 
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4.  Summary  and  Conclusions 


Impact  orientation  angle  histograms  were  generated  for  SRMU  segments  of  four  different 
lengths,  using  realistic  6D0F  free-fall  simulations.  For  high  initial  rotation  rates,  the  impact 
orientation  angle  was  found  to  be  uniformly  distributed  for  all  segment  lengths.  For  lower 
initial  rotation  rates,  longer  segments  were  found  to  trim  to  side-on  impact  in  slight  preference 
to  end-on.  For  the  shorter  segments  with  lower  initial  rotation  rates,  the  tendency  was  found 
to  be  toward  end-on  impact. 

A  method  for  computing  a  threshold  angle  for  side-on  impact  was  formulated.  For  impacts 
of  200  ft/sec,  this  angle  was  found  to  be  approximately  10°.  For  impacts  of  1000  ft/sec,  it 
was  found  to  be  approximately  23°.  For  impact  angles  lying  below  the  threshold,  the  yield 
for  the  event  is  taken  as  the  perfectly  side-on  value.  For  impact  angles  greater  than  twice  the 
threshold  value,  the  yield  for  the  event  is  taken  as  the  perfectly  end-on  value.  Intermediate 
cases  are  assigned  the  average  of  the  side-on  and  end-on  values. 
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Figure  1:  Lift  coefficient  for  16  ft  segment 


Figure  2:  Lift  coefficient  for  33  ft  segment 
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Figure  4:  Lift  coefficient  for  68  ft  segment. 
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Figure  10:  Moment  coefficient  for  33  ft  segment. 
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Figure  13:  Impact  angle  histogram  for  16  ft  segment  (apy  =  30  deg/ 
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Figure  14:  Impact  angle  histogram  for  16  ft  segment  ( —  3  deg / 


17 


0  10  20  30  40  50  60  70  80  90 

0  (deg) 

Figure  15:  Impact  angle  histogram  for  33  ft  segment  {apy  =  30  deg/sec) 
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Figure  16:  Impact  angle  histogram  for  33  ft  segment  ( =  3  deg/sec). 
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Figure  17:  Impact  angle  histogram  for  48  ft  segment  ( —  30  deg/sec) 
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Figure  18:  Impact  angle  histogram  for  48  ft  segment  (apy  =  3  deg/sec). 
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Figure  19:  Impact  angle  histogram  for  68  ft  segment  (apy  =  30  deg/sec). 
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Figure  20:  Impact  angle  histogram  for  68  ft  segment  (apy  =  3  deg/sec). 
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Figure  21:  Impact  orientation  of  cylinder. 
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